Proc. Indian Acad. Sci. (Math. Sci.) Vol. 113, No. 3, August 2003, pp. 321-332. 
Printed in India 



Probabilistic representations of solutions to the heat equation 

B RAJEEV and S THANGAVELU 

Indian Statistical Institute, R.V. College Post, Bangalore 560 059, India 
E-mail: brajeev@isibang.ac.in; veluma@isibang.ac.in 

MS received 16 November 2002 

Abstract. In this paper we provide a new (probabilistic) proof of a classical result 
in partial differential equations, viz. if is a tempered distribution, then the solution of 
the heat equation for the Laplacian, with initial condition <p, is given by the convolution 
of <p with the heat kernel (Gaussian density). Our results also extend the probabilistic 
representation of solutions of the heat equation to initial conditions that are arbitrary 
tempered distributions. 
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1. Introduction 

Let (X,),>o be a (i-dimensional Brownian motion, withXo = 0. Let <p G J?' (Mr), the space 
of tempered distributions. Let <p t represent the unique solution to the heat equation with 
initial value <p, viz. 

d,(p, = ^A(p, < t < T\ (po = (p- 

It is well-known that <p r = <p * p t , where p t {x) — ^^ rf / 2 e~(l*l 2 / 2f ) and V denotes convo- 
lution. When <p is smooth, say <p £ 5? ', the space of rapidly decreasing smooth functions, 
then the probabilisitc representation of the solution is given by the equality (p(t,x) = 
E(p(X t + x) and is obtained by taking expectations in the Ito formula 

<p (X, +x) = cp (x) + [ V(p (X s + x) ■ dX s + ^ / A(p {X s + x)ds. 
Jo 2 Jo 

Such representations are well-known (see 1 11213 141 ) and extend to a large class of initial 
value problems, with the Laplacian A replaced by a suitable (elliptic) differential operator 
L and (X t ) being replaced by the diffusion generated by L. A basic problem here is to 
extend the representation to situations where <p is not smooth. 

The main contribution of this paper is to give a probabilistic representation of solutions 
to the initial value problem for the Laplacian with an arbitrary initial value <p £ 5?' . This 
representation follows from the Ito formula developed in |9|, for the ^'-valued process 
(ix, <p), where T x (p is the translation of <p by x £ W 1 . Our representation (Theorem 2.4) 
then reads, % = Etx,<P where of course <p, is the solution of the initial value problem for 
the Laplacian, with initial value <p e S". In particular, the fundamental solution p t (x — •) 
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has the representation, Pt(x — ■) = Exx,8 x . However, the results of j9| only show that if 
<p G S?"p, then there exists q > p such that the process (xx, <p) takes values in y' q . Here for 
each real p, the y p s are the 'Sobolev spaces' associated with the spectral decomposition 
of the operator \x\ 2 — A or equivalently they are the Hilbert spaces defining the countable 
Hilbertian structure of y' (see (6)). 5?L, the dual of y p , is the same as y~ p . Clearly 
it would be desirable to have the process (xx,(p) take values in y' whenever (p G y'„. 
Such a result also has implications for the semi-martingale structure of the process (xx,) 
- it is a semi-martingale in y' j (Corollary 2.2) and fails to have this property in y' q for 
q < p + 1 (see Remark 5.2 of [5 1). 

Given the above remarks and the results of |9|, the properties of the translation oper- 
ators become significant. We show in Theorem 2.1 that the operators x x : y p — > y p for 
x G Mr, are indeed bounded operators, for any real p, with the operator norms being 
bounded above by a polynomial in |jc|. The proof uses interpolation techniques well- 
known to analysts. Theorem 2.4 then gives a comprehensive treatment of the initial value 
problem for the Laplacian from a probabilistic point of view. 

2. Statements of the main results 

Let (£1,^, (^t)t>0iP) be a filtered probability space with a filtration (J^,) satisfying 
usual conditions: & t = [\ s>t ^ s and J^b contains all P-null sets. Let (X t ) t >o be a d- 
dimensional, (^)-Brownian motion with Xo = 0. 

.y denotes the space of rapidly decreasing smooth functions on M. d (real valued) and 
<y' its dual, the space of tempered distributions. We refer to II II for formal definitions. 
For x G M d ,5t G <5f' will denote the Dirac distribution at x. Let {x x : x S M. d } denote 
the translation operators defined on functions by the formula x x f(y) = f(y — x) and let 
Xx : y —> y' act on distributions by 

(X x cp,f) = ((p,X- x f). 

The nuclear space structure of J?' is given by the family of Hilbert spaces ^ p ,p G K, 
obtained as the completion of under the Hilbertian norms {|| • || p } p eK defined by 

\\<p\\ 2 p ^(2\k\+d) 2p {<pM)\ 

k 

where <p G S?, and the sum is taken over k = {k\,... ,k c i) G Z+, \k\ = (k\ H hfed),{<p,%) 

denotes the inner product in L 2 (R d ) and {h k ,k G Z d + } is the ONB in L 2 (R d ), constructed 
as follows: for x = (xi,... ) Xd),hi c (x) — (x\) . . .h^Axa). The one-dimensional Hermite 

functions are given by ht{s) — Tj^ityTjz ^^' ''^Hi(s), where H((s) = (— l)V' 2 ^e~ 4 ' 2 
are the Hermite polynomials. While we mainly deal with real valued functions, at times 
we need to use complex valued functions. In such cases, the spaces <5^ p are defined in a 
similar fashion as above, i.e. as the completion of 5? with respect to || ■ || p . However, in the 
definition of ||<p||p above we need to replace the real L 2 inner product (<p,/ty) by the one 
for complex valued functions, viz. («p, yr) = f R d (p(x)\j/(x)dx and (cp,^) 2 is replaced by 
\((p,h k )\ 2 . It is well-known (see |6|7)) that y = f) p y p ,y" = \J p y p and y' p =: dual of 
y p — y^p. We will denote by (•,•)/>> the inner product corresponding to the norm || • \\ p . 

Let (Y t )t>o be an ^-valued, locally bounded, previsible process, for some p G K. Let 
dt : y p — > y p -i/2 be the partial derivatives, 1 < i < d, in the sense of distributions. Then 
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since <9,-, 1 < i < d are bounded linear operators it follows that (diY t ) t >o is an S^ p _\n- 
valued, locally bounded, previsible process. From the theory of stochastic integration in 
Hilbert spaces 1 8 1, it follows that the processes 

V s dX.A ,( fdiYsdX! 



r>0 \ JU / t>0 

are continuous j£> local martingales for 1 < i < d, with values in 5^ p and ^ p _i/2 respec- 
tively. If X t — (X, l ,...,Xf) is a continuous M^'-valued, J^ r -semi-martingale, it follows 
from the general theory that the above processes too are continuous ^-semi-martingales 
with values in and S? p _ \ h respectively. 

Theorem 2.1. Let pgR. There exists a polynomial P^-) of degree k — 2([|/?|] + 1) such 
that the following holds: For x £ R d , T x : ,5^ p — > ,5^ p is a bounded linear map and we have 

\\%Mp<Pki\A)\\<?\\v 

for all (p £ y p - 

In (|9|, Theorem 2.3) we showed that if (X t ) t >o is a continuous, rf-dimensional, & t - 
semi-martingale and <p £ S^ p C 5?' , then the process (Tx t <p)t>o is an ^-valued continu- 
ous semi-martingale for some q < p. Corollary 2.2 below says that we can take q = p — 1 . 

COROLLARY 2.2. 

Let (X t ) t >Q be a continuous d-dimensional, ^ t -semi-martingale. Let (p £ y p ,p £ ffi. Then 
(Tx,<p)t>o ^ an y p -valued, continuous adapted process. Moreover it is an 5f p -\-valued, 
continuous ^ t -semi-martingale and the following Ito formula holds in ,5f p - \: a.s. , V t > 0, 

d ,- t 

lx< (p = Tx (p - £ / di{xx s (p)AXi 
i=l J() 

\ t fdfcwWA (2.1) 



2 u 1 

where X t = (X t l , . . . ,Xf) and ((X 1 ,Xi) t ) is the quadratic variation process between (Xj) 
and (X/), 1 < i,j < d. 

Proof. From Theorem 2.1, it follows that (Tx,<p) is an ^-valued continuous adapted 
process. By Theorem 2.3 of |9), 3 q < p, such that (tx,<p) is an S" q semi-martingale and 
the above equation holds in S?q. Clearly each of the terms in the above equation is in 
y p -\ and the result follows. D 

The next corollary pertains to the case when (X t ) = {X} , . . . ,X t d ) is a cZ-dimensional 
Brownian motion, Xq = 0. In (|5|, Definition 3.1), we introduced the notion of an =5^p(= 
y-p,p > 0)-valued strong solution of the SDE 

&Y,= l -A{Y t )&t + VY t -dX tl 

Yo = <p, (2.2) 

where V = (d { , . . . , d d ) and A = £f =1 df. There we showed that if <p £ then the above 
equation has a unique J^'-valued strong solution, q > p + 2. Theorem 2. 1 implies that we 
indeed have an (unique) J^'-valued strong solution. 
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COROLLARY 2.3. 

Let <p G S^L. Then, eq. (2.2) has a unique 5^" p -valued strong solution onO <t <T. 

Proof. By Corollary 2.2, the process (Tx,<p), where (X t ) is a li-dimensional Brownian 
motion, Xq = 0, satisfies eq. (2.1). Further, 

Uniqueness follows as in Theorem 3.3 of 0. □ 

We now consider the heat equation for the Laplacian with initial condition <p S S^ p , for 
some p € K. 

d,(p, = ^A(p, 0<t<T, 

<Po = <P- (2.3) 

By an ^-valued solution of (2.3), we mean a continuous map t — > % : [0, T] — > 5 P such 
that the following equation holds in S?„- \ : 



<Pt = 9+ xA^di. (2.4) 

Jo 2 

Let {/i^ 1 } be the ONB in given by /if -1 = (2|Jt| +d)-^-^h k . We then have for 

;?<0andf<r: 

1*1=0 



1*1=0 



£ (<p ! /,r i )2_ 1 +2/ (9, J *jpv 1 d<9„Arv ] 



1*1=0 ■ /0 \ z 



ds 



1 

o \2 Z 



l<Pllp-i+2 /. (^A(p s ,(p s ) ds. 



p-i 

It follows from the results of [5 1 (the monotonicity condition) that for p < 0, 



1 



2 



2(-A<p,<p) +£||a i< p||2_ 1 <c 11911^ 
\ z /p-1 1=1 

for some constant C > for all 9 £ J?p. We then get 
\\<Ptf P -i <\\<pf P -i + C JjWsfp-ids. 

Hence for the case p < 0, uniqueness follows from the Gronwall lemma. Uniqueness for 
the case p > 0, follows from uniqueness for the case p < and the inclusion S^ p C Sf q for 
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q < p. It is well-known that the solutions of the initial value problem (2.3) in S^'(M. d ) are 
given by convolution of <p and p t (x), the heat kernel. That these coincide (as they should) 
with the ^-valued solutions follows from the 'probabilistic representation' given by 
Theorem 2.4 below. Define the Brownian semi-group (7)) r >o on 5? in the usual manner: 

T t <p(x)=<p*p t (x) f>0, T <p = <p 

where Pt(x) — p^57? e ^*^ 2 ^i* > an d '*' denotes convolution: f*g{x) — f R d f(y)g(x— 
y)dy. In the next theorem we consider standard Brownian motion (X t ). 

Theorem 2.4. (a) Let <p e ,S^ P . Then for t > 0, the S^p-valued random variable Tx,<p is 
Bochner integrable and we have 

E T X ,(P = <p*p t = T,(p. 

In particular, for every pSl, and T > 0,sup r<r \\T t \\ < °° where \\T t \\ is the operator 
norm ofT t : — > JJ'p. 

(b) For (j? S o5^ p , ?/ze initial value problem (2.3) /zas a unique S^p-valued solution <p, g/ven 
fry 

9, =Et Xi (P- 
Further <p, — > <p strongly in S^ p as t — > 0. 

3. Proofs of Theorems 2.1 and 2.4 

The spaces can be described in terms of the spectral properties of the operator // 
defined as follows: 

/// = (|x| 2 -A)/, /e^. 

If {/ijj is the ONB in L 2 (M. d ) consisting of Hermite functions (defined in §2), then it is 
well-known (see 1 10 1) that 

Hh k =(2\k\+d)h k . 

For / € 5?, define the operator H p as follows: 

HPf = Y,(2\k\+d) p (f,h k )h k . 

k 

Here p is any real number. For / G 5? and z = x+ iy € C define ff z / = La(2|^| + 
d) z (f,h k )h k and note that, ff z / = H^H^'f) = H^{H x f) and ff^ :L 2 ^L 2 is an isometry. 
Further, 

\m\\ 2 Q =£(2\k\+d) 2 *(f,h k ) 2 

k 



The following propositions (3.1, 3.2 and 3.3) may be well-known. We include the proofs 
for completeness. 
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PROPOSITION 3.1. 

For any p andq, \\H p (p\\ q - p = \\(p\\ q for (p G 5? . Consequently, H p : 5? q —* Yq-p extends 
as a linear isometry. Moreover, this isometry is onto. 

Proof. Let h\ = (2\k\+d)~ p h k . Then from the relation ((p,h k ) p = (2\k\+d) 2p {(p,h k ) it 
follows that {h p k } is an ONB for y p . Let (p e.Y. Since 

HP(p = Y J (<P,h k )(2\k\+dyh k 

k 

= ^{ ( p,h k ){2\k\+d)' , hi- p , 

k 

weget||^(p||2_ p =||<p||2. 

To show that H p is onto, consider y/ G y q - p , 

V = L(Y,hr P ) q -ph q k - p . 

k 

Defining <p =: Y.k(V,h k ~ p ) q - P hl, we see that <j? G y q . Also, 

k k 

Let Aj=Xj + dj and A+ = Xj — dj, 1 < j <d. Then it is easy to see that 

Z 7=1 

For multi-indices a = (0Ci, . . . , a^) and J3 = (/Ji, . . . ,/3 d ) we define 

A°=:A?...A?, (A+)^:(A+)ft...(A+)fe. 
For an integer £ > and x G M, recall that 

h e (x) = - - ] xl/0 e-< x2 /2) tf ( x ) 

where is the Hermite polynomial defined by 

= (-ljV^e-* 2 . 
It is easily verified that 

(e-^(x)) =«(e-^_i(x)) , 
'x-±) (e-^(x))=e"^ +1 (x). 



It then follows that 



A+fyfa) = yj2(kj + l)h kj+1 (Xj), 

Ajh kj (xj) = y/2kjh k .-i(xj). 
Iterating these two formulas we get the following: 
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A a h k (x)=2^ 2 



k\ 



(k-a)l 



1/2 



where k\ = k\\ . ..kjl. 
PROPOSITION 3.3. 

For all m > 0, 3 constants C\ — C\ (m) and C2 = Ci(m) such that the following hold: 

(a) For all fey, 

||/||™<Ci £ ||A«(A + )^/|| <C 2 ||/|| m . 
|a|+|j3|<2m 

(b) For all fey, 



||/||»<Ci £ lk a ^/||o<c 2 ||/|U. 

|a|+|j3|<2m 
Proo/ (a) We can write 

# m = E C ali A a (A+)P, 

\a\ + \p\<2m 

where C a p are constants. Since ||/||m = ||/? m /||o, the first part of the inequality follows. 
To show the second half of the inequality it is sufficient to show that for / S and 

\a\ + \p\ < 2m, \\A a (A+)PH- m f \\o < C a/3 1|/|| . Now, 



\\ A «{A+fH- m fg = Y,{A a {A + fH- m f,h t ) 2 



-1 2 



l^(2\k\+d)- m {f,h k ){A a (A + fh k M) 



Z(2\k\+d)- m (f,h k }C^ !a (h k+ p_ a ,h f ) 



Y,(2\£ + a-l5\+dr 2m Cl a _^^ a (f,h e+a _p) 2 , 



where the sum is taken over £ = {l\ , . . . , ij) such that lj + ccj — /?,■ > for 1 < j < d and 
where we have used Proposition 3.2 in the last but one equality above. From the same 
proposition, it follows that 

(2\a + e-i3\+d)- 2m C 2 +a _ p ^ !a 



are uniformly bounded in I for \a\ + |j3 1 < 2m and the second inequality in (a) follows. 
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(b) Since = ||//"'/||o and clearly H m = L|a|+|/3|<2mC a /3Jc a <3' 3 , the first inequality 
follows. To prove the second inequality, note that 



1,. , 1 



Xj = ^(Aj+Af), dj = -(Aj-Aj). 
Hence, using [Aj,A^] = 5,jJ, 

x a dP = £ CuA k (A + ) e 

\k\+\i\<\a\+\p\ 

and hence by part (a) we get 

£ ||*°^/||o<C, £ ||A*(A+)V||o<C2||H»/|| . 

|a| + |0|<2m \k\+\e\<2m 

□ 

Proof of Theorem 2.1. We first show that for an integer m > 0, 

hxf\\m <P2m(\A 



where P2m{t) is a polynomial in t £ M. of degree 2m with non-negative coefficients. This 
follows from Proposition 3.3: 

IMU<C, £ \\y a d^ x f\\ Q 

\a\+\p\<2m 

<d £ ||(y + x) a ^/|| . 

a|+|/3|<2m 

The last sum is clearly dominated by ^md^DII/IL f° r some polynomial i^m- If m < p < 
m + 1, where m > is an integer, we prove the result using the 3-line lemma: for /, g G & 
let 

F(z) = (Hh x H-'f,g) Q . 

Then from the expansion in L 2 for the RHS it is verified that F(z) is analytic in m < 
Re z < m + 1 and continuous in w < Re z < '« + 1 . We will show that 

|F(m + o,)|<P 2ffl (|x|)||/|| ||g||o, 

\F(m+l+iy)\<P 2{m+1) (\x\)\\f\\ \\g\\o (3.1) 
for — °° < y < oo. Hence from the 3-line lemma 1121 . it follows that 

l^(p+v)l<(^(W)ll/l|olkllo)^ 1 - p (^i ) (W)ll/llolkllo)^ m 
<flfc(W)ll/llolkllo, 

where Pk{t) is a polynomial in t of degree k — 2([p] + 1). It follows that 

IK/|| P </KW)II/IU- 
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Using the fact that J^_ p = ,9" p we get \\T x f\\- p < Pk(\x\)\\f\\- P for m < p < m+ 1. 
The following chain of inequalities establish the inequalities (3.1): 

\F(m + ty)\ < WH'^y^H-^fWoWgh 
< \\H m r x H-^f \\o\\g\\o 
= \\T x H- (m+,y) f\\ m \\ 8 \\o 

<P2,n(\ X \)\\H-^f\\ m \\g\\o 

= P 2m (N)||//-*/||o|k||o 

= ft™(W)ll/l|o|k||o. 



This completes the proof of Theorem 2.1. D 
Proof of Theorem 2.4. (a) Let (p £ y p ,p £R. From Theorem 2.1 we have 

hxM\p<Pk(\Xt\)\\<p\\p, 



where P% is a polynomial. Since EP^{\X t \) < Bochner integrability follows. For \\f £ 



y/,J T x (pp t (x)dxj = J {\i/,T x (p)p t (x)dx 

= J p,(x)6x J y{y)<p{y-x)dy 
= jy(y)&y J (p(y-x)p t (x)6x 
v(y)<p*pt(y)dy 



= (W,<P*Pt)- 

The result for <p £ follows by a continuity argument: Let (p n £ <p n — <p in S? p . 
Hence (p n * p t — > (p * p t weakly in 3" . Hence, 

(W,<P*Pt) = lim (W,%*Pt) 

= lim / \i/(y)(p n *p,(y)dy 



= lim / (\j/,T x (p n )p t (x)dx 
= J (y,i; x (p)p t (x)dx 
V, J t x <p pt(x)6xj , 



where we have used DCT in the last but one equality. That T t : ,5P p — ► Sf p is a (uniformly) 
bounded operator follows: 
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\T,<p\\p = \\<p*p,\\p = \\ET Xl <p\\p 



J T x <p p,(x)dx 



< / \\T x (p\\pp t (x)dx 



<\\(p\\ P J P k (\x\)p t (x)dx<C\\cp\\ p , 
where C = sup s<r J P/ ( (\x\)p s (x)dx < °°. 

(b) Let (X t ) be the standard Brownian motion so that (X',X^) = for i ^ j. Equation (2.1) 
then reads, for <p 6 ^ p ,p £ R, 



(3.2) 



T Xl (p = <p- V(Tx,<p)-dX 4 .+ - / A(T X ,(p)ds. 
J0 2 Jo 

The stochastic integral is a martingale in 5?p—\\ 



diizx^dX's 



= Q I J ||a ! (T^)||;_ 1 /?. v (x)dx)di 



<c 3 \\<p\ 



P k (\x\)p s {x)dxj ds 



Let <p ( = ETx,<P- Taking expected values in (3.2) we get eq. (2.4). Hence (p t is the solution 
to the heat equation with initial value <p € S^ p . The uniqueness of the solution is well- 
known and also follows from the remarks preceeding the statement of Theorem 2.4. 

To complete the proof of the theorem, we need to show that q>, — > <p in S^ p as t { 0. 
Let ^ denote the Fourier transform, i.e. .^f{£, ) = J e~'( x '^ f(x)dx for f £ ,Y. Then 
J^" extends to by duality, where we consider 5?' as a complex vector space. Since 
,^(h n ) = (—y/—l)"h„ (1 10 1, p. 5, Lemma 1.1.3), & acts as a bounded operator from S^p 
to for all p. Let <p e J^ p . 

(p t -(p = T t (p-(p = (S t (&(p)), 

where 

S f <p(x) = &{T t -I)&- l <p{x) = ( e -('/ 2 )W 2 - l)<p(jc). 
Clearly, 5 f : ^ p — ► ^ is a bounded operator and 

||<ft-(p|| p = ||s f (jr<p)|| p . 

The following proposition completes the proof of the theorem. 
PROPOSITION 3.4. 

Let <p e y p ,p E K. TTzen ||5 f <p|| p — > as t — > 0. 
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Proof. We prove the proposition by showing that (i) S t : 5^ p — > ^ are uniformly bounded, 
< t < T and (ii) ||S r <p|| p — > for every <p e -V, as t — > 0. Let us assume these results for 
a moment and complete the proof. 

Let £ > be given. By (i), there is a constant C > such that 

sup \\S t f\\ p <C\\f\\ p ,feS' p . 

0<t<T 

Choose <p € ^, so that ||/- <p|| p < (^). Then, 

\\S t f\\p<\W-<p)\\ P +\\S t <p\\p 
<e/2 + \\S t (p\\ p . 

Now choose 5 > such that \\S t (p\\ p < e/2 for all < f < 5, to get ||5 r /|| p < £ for all 
0<r<5. 

Since S, = &{T t - 1) J?" 1 , (i) follows from the fact that 7} : S? p -> ^ p are uniformly 
bounded (Theorem 2.4a) and ^ : ^ p — » =5^ is a unitary operator. The proof of (ii) is by 
a direct calculation when p = m is a non-negative integer. 

\\S t ( P \\m = \\H m S t ( P \\o<C 1 £ Hjc^S^Ho. 

|a|+|j3|<2m 

Since S t <p{x) = (e^'/ 2 ^ 2 - l)<p(x), by Leibniz rule 

||* a <?^<p||o< £ C Mr ||* a ^(e-™ 2 -l)d>|| . 

iMl+lrhlPI 

When 7^ 0, we have 

|| x a^ (e -('/2)W 2 _ 1)5 r (p || < C2f lMl|| (p || m 

and when n = 0, using the elementary inequality 1 1 — e~"| < C3M, u > we get 

IK(e- (t/2)W2 -l)a>|| <C 4 ?||(p|| m+ i. 

Therefore, ||S,<p|| m < C t ||<p|| m +i for some constant C, which shows that ||S,<p|| m — > as 
f — > 0. If /? is real and m is a non-negative integer such that p <m,we have 

||5 f (p|| P <||5 ( (p|| ffl <Cf||(p|| m+1 

and so ||S f <p|| p — > as f — > in this case as well. □ 
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